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PREFACE 


This  report  was  prepared  by  John  B.  Howe  of  the  Mount  Washington  Observatory,  Gorham,  New 
Hampshire.  Report  preparation  was  funded  through  DA  Project  4A161  \02KT2A,  Re  search  in  Snow, 
Ice  and  Frozen  Ground,  Task  FS,  Work  Unit  005,  Concepts  of  Mesoscale  Winter  Meteorological 
Processes:  some  of  the  work  described  was  funded  through  other  CRREL  projects  during  the  past  1 5 
years.  Other  studies  by  the  author  discussed  in  this  report  were  funded  by  USAF  Air  Materiel 
Command.  Sigma  Xi,  and  the  National  Science  Foundation. 

Many  past  and  present  staff  members  of  the  Mount  Washington  Observatory  and  of  the  Aeronau¬ 
tical  Icing  Research  Laboratories  have  been  involved  in  gathering  the  data  for  the  many  hundreds  of 
Rotating  Multicylinder  exposures  that  the  author  has  analyzed  or  examined.  Otherindividuals  whose 
counsel  has  been  influential  include  R.  Cunningham,  formerly  of  Air  Force  Cambridge  Research 
Laboratories:  W.  Howell,  formerly  of  the  Mount  Washington  Observatory;  P.  Perkins,  formerly  of  the 
National  Advisory  Committee  for  Aeronautics/National  Aeronautics  and  Space  Administration;  and 
C.  Franklin,  formerly  of  the  Aeronautical  Icing  Research  Laboratories.  The  support  and  advice  of  S. 
Ackley,  J.  Govoni,  K.  Jones,  C.  Ryerson.  and  W.  Tucker  of  CRREL  are  gratefully  acknowledged. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional  purposes.  Citation  of 
brand  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial 
products. 
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NOMENCLATURE 


a  droplet  radius  (cm) 

C  cylinder  diameter  (cm) 

D  droplet  diameter  (urn) 

E  collection  efficiency,  dimensionless 

£m  overall  weighted  collection  efficiency 

f\  9(V'/n)  (pa)2 


K 
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LWC 

^Re 
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vv 

vr 
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<t> 

Subscripts 

a 

av 

d 

i 

o 

w 


pwDV  ..  .  , 

inertia  parameter,  — - ,  dimensionless 

9ri  C 

length  of  ice  collection  (cm) 

liquid-water  content  (g/m3)  (informal  usage) 

free-stream  Reynolds  number  with  respect  to  droplet 

diameter,  paDV7ri,  equivalent  to  (/(0)l/2 

exposure  duration  (s) 

free-stream  airspeed  (cm/s) 

liquid-water  content  (g/m3) 

cylinder  weight  (g) 

absolute  viscosity  of  air  (g/cm  s) 

density  (g/cm3) 

(Wrc)2  _9p \CV 


K 


PwH 


,  dimensionless 


air 

average 
volume  median 
with  ice  collection 
without  ice  collection 
liquid  water 


Rotating  Multicylinder  Method  for  the  Measurement  of 
Cloud  Liquid- Water  Content  and  Droplet  Size 

JOHN  B.  HOWE 


INTRODUCTION 

The  Rotating  Multicylinder(RMC)  method  has  been 
in  use  since  the  late  1940s.  but  no  description  of  it  has 
ever  appeared  in  the  open  literature,  although  it  is 
among  the  best,  and  is  certainly  the  simplest  and  most 
reliable,  of  any  technique  for  the  measurement  of  cloud 
liquid-water  content  and  droplet  size  in  icing  condi¬ 
tions.  This  report  attempts  to  provide  sufficient  infor¬ 
mation  to  enable  someone  unfamiliar  with  the  method 
to  build  and  use  the  instrument.  Included  also  is  a  review 
of  the  theory  upon  which  the  data  analysis  procedure  is 
based. 


HISTORY 

The  RMC  method  was  developed  at  the  Mount 
Washington  Observatory  (MWO)  between  1940  and 
1945  after  prior  use  of  a  single  rotating  collector  for 
measuring  liquid  water  content  (LWC)  of  clouds 
(Arenberg  1939,  1941).  The  method  was  given  a  rigor¬ 
ous  theoretical  basis  by  Langmuir  ( 1 944)  and  Langmuir 
and  Blodgett  1946).  Langmuir  ( 1944)  mentions  some 
prior  work  by  F.  Albrecht  in  Germany  and  by  M. 
Glauert  in  England.  Other  workers  involved  in  the 
development  of  the  method  at  the  MWO  included  V. 
Clark  of  the  MWO  staff.  R. Cunningham  of  MIT,  and  V. 
Schaeferof  General  Electric.  During  the  late  1940s  and 
the  1950s  the  method  was  used  extensively  for  icing 
studies  in  natural  wind,  wind  tunnels,  and  aircraft  by 
various  groups,  including  the  National  Advisory  Com¬ 
mittee  for  Aeronautics  ( N  AC  A),  the  Aeronautical  Icing 
Research  Laboratories  tAIRL),  and  the  MWO.  Tabu¬ 
lation';  of  data  collected  with  the  method  appeared  in  a 
series  of  quarterly  reports  by  the  MWO  staff  ( 1 949-5 1 ). 
Ambrosio  (1950)  tabulated  all  results  obtained  by  all 
agencies  during  the  period  1 945  1950,  a  total  of  some 
3600  observations. 


METHOD 
General  description 

A  cylinder  exposed  for  a  known  period  of  time  in  an 
airstream,  with  its  axis  perpendicular  to  the  flow,  will 
“sweep  out”  a  known  volume  of  air.  If  the  air  contains 
droplets  of  supercooled  water,  it  might  be  assumed  that 
all  the  droplets  in  the  swept-out  volume  would  strike 
and  freeze  on  the  cylinder.  This  assumption  is  nearly 
correct  when  the  droplets  are  very  large  and/cr  the  air¬ 
speed  is  very  high.  However,  as  the  droplet  size  or  air¬ 
speed  becomes  smaller  the  droplets  have  less  momen¬ 
tum,  and  those  near  the  edges  of  the  swept-out  volume 
will  tend  to  follow  the  airstream  lines  around  the  cylin¬ 
der  without  striking  it.  If  the  diameter  of  the  cylinder  is 
increased,  this  effect  becomes  more  pronounced.  The 
ratio  of  the  mass  of  liquid  water  striking  the  cylinder  to 
the  total  mass  of  liquid  water  in  the  swept-out  volume  is 
called  the  collection  efficiency  E  of  the  cylinder.  Col¬ 
lection  efficiency  varies  directly  with  airspeed  and 
droplet  size  and  inversely  with  cylinder  diameter.  Col¬ 
lection  efficiency  is  affected  to  a  smaller  degree  by 
changes  in  the  density  and  viscosity  of  the  air. 

The  liquid-water  content  ofacloud  (LWC  forinfor- 
mal  usage,  W  in  mathematical  usage)  is  the  mass  of 
liquid  water  per  unit  volume.  The  mass  of  ice  collected 
on  acylinder  per  unit  of  swept-out  volume  is  equal  to  the 
LWC  of  the  cloud  multiplied  by  the  collection  efficien¬ 
cy  of  the  cylinder,  or  EW.  In  practice.  EW  for  a  cylinder 
is  calculated  from  the  relation 

weight  of  collected  ice 

fcvv  -  - — — - — - 

swept-out  volume 

where  the  swept-out  volume  is  calculated  from  mea¬ 
sured  values  of  airspeed,  exposure  duration,  cylinder 
length,  and  average  cylinder  diameter.  (Average  diam¬ 
eter  is  the  mean  of  the  diameters  at  the  beginning  and 
end  of  the  exposure,  the  cylinder  having  been  rotated 
during  exposure  to  produce  a  cylindrical  ice  collection.) 


In  the  RMC  method,  several  cylinders  of  different  diameter  of  the  cylinders  is  not  critical  and  standard 

diameters  are  exposed  simultaneously;  the  resulting  ice  aluminum  tubing  sizes  maybe  used  such  that  each  cyl- 

collections  are  measured  and  the  values  of  cylinder  mder  is  approximately  twice  the  diameter  of  its  smaller 

diameter  vs  EW  are  plotted  on  logarithmic  paper.  neighbor.  A  slightly  loose  fit  between  all  components 

Theoretical  studies  of  droplet  trajectories  nearcylin-  makes  for  ease  of  disassembly, 
ders  have  resulted  in  graphs  of  a  dimensionless  param-  The  bottom  endof  the  shaft  is  threaded  to  fit  the  drive 

eter  1  IK  vs  cylinder  collection  efficiency;  the  expert-  motor  shaft.  At  MWO  we  use  gear  motors  giving  about 
mental  data  plot  is  fitted  to  one  of  these  curves,  LWC  is  1  or  2  rpm;  rotation  speed  should  be  close  to  this  range, 

read  off  directly,  and  a  simple  calculation  yields  the  but  is  not  critical.  (Aircraft  speeds  require  faster  rota- 

droplet  diameter.  tion.)  We  have  several  taper  mounts  on  our  tower  railing 

so  that  the  drive  motor  with  the  RMC  can  be  exposed  on 
Instrumentation  whatever  happens  to  be  the  windward  side.  Provision 

Figure  1  shows  an  instrument  assembled  and  disas-  may  have  to  be  made  for  adjusting  the  angle  of  the  motor 

sembled;  working  drawings  and  a  detailed  description  mount  so  that  the  instrument  shaft  is  perpendicular  to 

are  given  in  Appendix  A.  Manv  other  designs  have  been  the  wind  vector.  Before  any  exposures  are  made,  the  site 

used,  but  this  one  is  the  simplest  for  natural  wind.  The  should  be  checked  with  a  straight  rod  collector  to  assure 
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a  Assembled. 

Figure  I .  Rotating  nmlticylinder  instrument. 


I).  Disassembled. 
Figure  1  (cant'd) 


that  the  instrument  will  be  outside  the  boundary  layer 
with  uniform  impingement  along  the  collector. 

Wind  speed  must  be  measured  as  accurately  as  pos¬ 
sible.  We  are  fortunate  in  that  the  MWO  official  ane¬ 
mometer  is  a  heated,  vane-mounted  Pitot-static  tube.  If 
an  anti-iced  anemometer  is  not  available,  an  accurate 
cup  anemometer  can  be  used  if  it  is  put  in  place  just  be¬ 
fore  the  beginning  of  an  RMC  run.  since  the  small 
amount  of  ice  or  rime  that  collects  on  the  cups  during  the 
exposure  should  not  be  detrimental. 

Weighing  (to  ±0.01  g).  length  measurement  (to 
±0.05  cm),  and  timing  (normally  to  ±5  s)  require  no 
special  equipment.  Ice  diameter  is  most  easily  mea¬ 
sured  with  a  inicrometercaliper.  Care  must  be  taken  not 
to  depress  the  rime  surface  when  measuring  a  diameter: 
the  caliper  must  be  adjusted  by  eye.  net  by  feel.  For  the 
three  largest  cylinders  this  measurement  is  most  easily 
done  by  measuring  one  side  of  the  iced  cylinder  over  an 
anvil.  (A  device  for  this  is  easily  made  by  inserting  a 
length  of  2-mm  wire  vertically  into  a  small  piece  of 
plywood.) 

Exposure  procedure 

A  step-by-step  outline  of  exposure  and  data-reduc- 
tion  procedure  is  given  in  Appendix  C. 

On  Mt.  Washington  we  make  exposures  in  natural 
winds  ranging  from  -30  to  - 1 50  km/hr.  For  the  highei 
speeds  of  aircraft  it  may  be  advisable  to  use  larger  cyl¬ 
inder  diameters;  for  ground  sites  w  ith  winds  of  less  than 
30  ktn/hr  it  may  be  necessary  to  employ  a  carefully 
designed  small  wind  tunnel. 


The  instrument  must  be  at  a  subfreezing  temperature 
prior  to  the  beginning  of  an  exposure.  Exposure  dura¬ 
tions  in  Mt.  Washington  conditions  range  from  5  to  40 
minutes.  Our  rule  of  thumb  is  to  end  an  exposure  w  hen 
the  collection  on  the  smallest  cylinder  has  reached 
approximately  the  original  diameter  of  the  no.  2  cy  lin¬ 
der.  This  is  a  compromise  between  the  desire  to  mini¬ 
mize  percentage  weighing  errors  and  the  need  to  main¬ 
tain  the  spread  in  cylinder  diameters.  Disassembly  and 
measurement  after  the  exposure  is  usually  a  straightfor¬ 
ward.  albeit  rather  delicate,  job.  but  when  hard  clear  ice 
is  collected  it  helps  to  have  a  warm  cutting  tool.  (A  light 
bulb  in  an  open-ended  tin  can  w  ith  sheet  copper  blades 
soldered  to  one  end  is  ideal.)  Disassembly  and  length 
and  diameter  measurement  must  be  done  in  a  coldroom. 
Weighing  will  usually  be  done  at  the  same  time,  but  if 
this  is  inconvenient  the  iced  cylinders  can  be  stored  in 
canisters  for  later  weighing.  If  part  of  the  collection  on 
a  cylinder  is  accidentally  broken  off.  the  remaining  ice 
can  be  trimmed  to  give  a  “square"  end:  the  length  di¬ 
mension  required  is.of  course,  that  of  the  ice.  not  that  of 
the  cylinder  itself. 

Data  reduction — preliminary 

Several  quantities  must  be  computed  prior  to  the 
curve-matching  step: 

1.  Average  diameter  (Cav.  cm)  for  each  cylinder,  the 
mean  of  the  iced  and  bare  diameters. 

2.  Product  of  LWC  andcollection  efficiency  <£31  ' 
for  each  cylinder. 


3 


(1) 


EW  =  v'-i  x  106,  e/nv’. 
VLTC  av 


3.  Air  density  (pa) 


p|  =  0:348  (jxessu,-efntb[)x  ,Q-3  g/cm-3  (2) 
(temperature  [K  |) 


4.  Airspeed  divided  by  viscosity  (V/r|) 

11  =  - __ L: - x  io6.  (3) 

n  2.48  (temp[K])° 754  S/cm  s 

5.  An  intermediate  factor  (f.) 

/,  =9(\7n)(pa)2.  (4) 

The  values  of  C  are  then  plotted  vs  ElTon  three-cycle 
logarithmic  tracing  paper  (Fig.  2). 

Data  reduction — curve  matching 

The  theoretical  curves,  to  which  we  fit  the  plots  of 
Cav  vs  EW.  are  graphs  of  a  dimensionless  parameter 
( VK)  vs  cylinder  collection  efficiency  (EM),  drawn  for 
selected  values  of  Kt).  (Kb  is  a  dimensionless  parameter 
equal  to  the  square  of  the  Reynolds  number  for  the 
volume-median  droplet  size.)  For  each  value  of  £<{>  there 
is  a  family  of  curves  labeled  "A”  through  “J"  ("1”  is 
omitted).  The  "A"  curves  assume  all  droplets  are  the 
same  size,  while  the  others  assume  progressively  broad¬ 
er  droplet  size  distributions.  The  curves  must  be  drawn 
on  the  same  type  of  graph  paper  used  for  plotting  the 
data  points.  It  is  convenient  to  put  all  nine  curves  of  each 
family  on  one  sheet  of  paper,  with  the  vertical  axis 
shifted  slightly  foreachcurve.  A  small  invertedT  marks 
the  point  \/K  =  I .  E  —  1  for  each  curve.  An  example  of 
these  curves  is  shown  and  the  theory  of  their  use  is 
described  and  explained  more  fully  in  Appendix  B. 

The  matching  of  the  C  vs  EWplot  to  the  theoretical 
curves  must  be  done  using  a  curve  whose  Kb  value  is 
closest  (on  a  logarithmic  scale)  to  the  actual  value  for 
the  exposure  being  analyzed:  this  value  is  of  course  not 
known  a  priori,  since  Kb  is  a  function  of  droplet  size. 
Therefore,  one  must  make  a  guess  and  try,  e.g..  the  Kb 
=  200  family.  (After  some  experience  one's  first  guess 
will  usually  turn  out  to  be  correct. )  The  procedure  is  as 
follows: 

1.  Slide  the  graph  paper  with  your  plotted  points 
around  on  the  theoretical  curves,  keeping  the  axes 
parallel,  until  you  find  the  best  match. 

2.  Read  on  your  plotted  graph  the  value  of  C  which 
lies  over  the  \/K=  I  axis  of  the  theoretical  curve:  this  is 
the  value  ofcylinderdiamcter  for  which  AT  =  1(<7K=1)  for 
the  conditions  of  this  exposure. 

3.  Multiply  this  number  times  the  value  oij\  (pre¬ 


viously  determined  by  eq  4)  to  find  the  value  of  AT0.  If 
this  is  closer,  on  a  log  scale,  to  200  than  it  is  to  2d  or  1 000 
you  have  used  the  correct  family  of  curves;  if  not,  repeat 
the  matching  process  using  the  appropriate  family. 

4.  Trace  the  curve  lightly,  and  trace  the  correspond¬ 
ing  inverted  T. 

5.  Now  read  the  value  ofLWC  for  your  exposure  on 
the  horizontal  scale  of  your  plot  over  the  £M  =  i  axis  of 
the  theoretical  curve  (in  line  with  the  inverted  T). 

6.  Use  the  value  of  Cav  forwhich  1  IK-  1  tocompute 
the  volume-median  droplet  diameter  (Dd)  and  Ko : 


10  ,  pm 


(5) 


Kb  — /)(C|£_| ).  (6) 

See  Figure  2  for  examples  of  curve  matching  for 
exposures  made  on  Mt.  Washington.  Figure  2a  is  from 
a  run  for  which  a  mistake  was  made  in  the  measurement 
of  the  no.  3  cylinder;  the  effect  of  the  mistake  is  obvious 
and  easily  allowed  for.  Figure  2b  is  typical  of  runs  made 
when  wind  speed  is  low  and/or  Dd  is  small;  there  was 
only  a  dusting  of  rime  on  the  tw'o  largest  cylinders, 
w  hich  were  not  measured,  and  the  no.  4  cylinder  had  a 
slightly  spotty  collection.  In  such  cases  points  for  cylin¬ 
ders  with  spotty  col’ections  are  best  ignored,  even 
though  it  may  be  possible  to  achieve  an  approximate  fit 
of  all  the  points  to  a  theoretical  curve  of  narrower  dis¬ 
tribution.  Figure  2c  illustrates  the  opposite  set  of  condi¬ 
tions,  with  high  wind  speed  and  large  Dd.  A  step-by-step 
explanation  of  the  exposure  and  analysis  procedure  is 
given  in  Appendix  C. 

Computer  processing 

All  the  computations  described  above  can  be  done 
by  a  very  simple  computer  program;  with  a  hand  calcu¬ 
lator  they  can  be  done  almost  as  easily.  Attempts  have 
been  made  to  computerize  the  curve-matching  process, 
but  in  my  opinion  this  is  a  mistake.  Aside  from  consid¬ 
erations  of  questionable  cost-effectiveness,  “wrong" 
matches,  due  to  measurement  errors  or  spotty  collec¬ 
tions,  may  result  and  will  be  difficult  or  impossible  to 
detect.  A  strong  recommendation  for  the  method  de¬ 
scribed  in  the  preceding  sections  is  that  such  errors  can 
usually  be  recognized  before  they  do  any  damage. 

Use  of  method  in  warm  fog 

The  RMC  can  be  used  in  above-freezing  tempera¬ 
tures  if  the  cylinders  and  intercylinder  adapters  are 
made  of  a  porous  material;  it  may  be  simpler  to  wrap  the 
regular  metal  parts  in  a  porous  substance  such  as  Milli- 
pore  filter  material.  Weights  must  be  measured  before 
as  well  as  after  each  exposure  and  great  care  is  needed 
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c.  typical  run  with  high  wind  speed  and  large  Dd  fW  =  0.68  glm3,  C(Kml)  =  33.0. 
Dj  =  39  urn). 

Figure  2  (cont'd).  Examples  of  curve  matching. 


to  ensure  that  the  cylinders  (particularly  the  smaller 
ones)  do  not  become  saturated. 


DISCUSSION  OF  ACCURACY 

This  discussion  will  begin  by  considering  the  accu¬ 
racy  of  an  “ideal"  exposure  resulting  in  complete  cov¬ 
erage  of  at  least  five  cylinders.  (Ideal  exposures  are  ac¬ 
tually  quite  common,  comprising  perhaps  half  of  all 
runs  made  at  MWO.)  Following  this  we  will  consider, 
first,  conditions  that  may  result  in  reduced  accuracy, 
and  then  discuss  hidden  problems  due  to  possible  errors 
in  the  theoretical  curves.  Finally  we  will  briefly  review 
the  results  of  comparison  tests,  both  between  RMCs  and 
between  the  RMC  and  other  instruments. 

“Ideal”  runs 

If  the  physical  measurements  have  been  made  with 
reasonable  care,  the  plotted  data  points  should  all  be 
within  5%  of  their  true  locutions.  (The  most  likely 
source  of  greater  error  is  in  the  measurement  of  wind 
speed.) 

An  ideal  exposure  will  result  in  an  almost  perfect 
curve  match;  if  the  curve  is  similar  to  that  of  Figure  2a 


accuracy  in  W  and  Dd  will  be  comparable  to  that  of  the 
original  measurements,  except  that  errors  in  wind  speed 
measurement  have  a  smaller  effect  on  the  computed  Dd 
(eq  5).  If  the  overall  slope  of  the  curve  is  low,  as  in 
Figure  2b,  accuracy  in  W  will  be  somewhat  lower  be¬ 
cause  there  is  less  constraint  on  the  “sideways”  position 
duri  ng  matching.  For  cases  like  Figure  2c,  with  a  steeply 
sloping  plot,  accuracy  in  Dd  will  be  lower  because  there 
is  less  constraint  on  the  “up-and-down"  position. 

Errors  due  to  carelessness  in  measurement  or  to 
breakage  of  ice  will  inevitably  occur,  but  will  rarely 
affect  more  than  one  or  two  cylinders.  Such  errors  can 
usually  be  easily  recognized  in  the  data  plot  and  allowed 
for  during  curve  matching,  so  accuracy  will  not  be 
seriously  affected  (Fig.  2a). 

Thus  for  ideal  runs,  assuming  the  validity  of  the 
theoretical  curves,  accuracy  in  W and  Dd  is  estimated  to 
be  ±10%  or  better.  This  includes  the  effect  of  small 
errors  caused  by  the  actual  value  of  K$  being  different 
from  the  value  used  to  construct  the  theoretical  curve. 

“Non-ideal”  runs 

When  Dd  is  very  small  and/or  wind  '■need  is  low,  the 
larger  cylinders  may  have  a  spotty  col  ction  (individ¬ 
ual  kernels  of  rime)  or  no  collection  a;  all.  Figure  2b, 
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referred  to  in  a  previous  section,  shows  an  example  of 
this.  Sometimes  there  will  be  measurable  collections  on 
only  three  cylinders,  rarely  on  only  two.  This  has  a 
potentially  drastic  effect  on  accuracy,  especially  since  a 
measurement  error  on  one  of  the  cylinders  cannot  be 
detected  in  the  data  plot. 

The  plotted  points  for  cylinders  with  spotty  collec¬ 
tions  tend  to  fall  below  the  proper  theoretical  curve. 
Though  they  can  often  be  forced  into  an  approximate  fit 
by  going  to  a  curve  for  a  narrower  distribution,  such  a 
forced  fit  will  be  of  questionable  validity.  A  notation 
should  always  be  made  on  the  measurement  work  sheet 
when  collections  are  spotty ,  even  when  coverage  is  near 
1 00%,  and  the  data  points  for  any  spotty  cylinders  must 
be  treated  with  suspicion  during  curve  matching.  In 
these  cases  the  reduction  in  accuracy  can  range  from 
insignificant  to  serious  and  will  have  to  be  judged  case 
by  case.  It  is  helpful  to  include  with  the  analysis  results 
for  each  run  a  notation  as  to  the  number  of  points  that  fit 
the  theoretical  curve  within  approximately  5%. 

Runoff 

Moderate  temperatures  and/or  high  rates  of  water 
impingement  can  result  in  some  fraction  of  the  collected 
water  running  off  without  freezing;  the  smaller  cylin¬ 
ders  are  most  likely  to  be  affected.  The  collection  will 
be  of  clear  ice,  although  the  occurrence  of  clear  ice  does 
not  necessarily  indicate  runoff.  If  possible,  the  instru¬ 
ment  should  be  checked  visually  during  exposure  when 
there  is  any  likelihood  of  this  phenomenon;  it  is  usually 
easy  to  see. 

Exposures  with  runoff  will  of  course  yield  only  a 
“lower  limit”  value  of  LWC  and  a  very  unreliable  value 
of  Dj,  even  if  accurate  measurement  of  the  (usually 
rapidly  melting)  ice  is  possible. 

Accuracy  of  theoretical  curves 

The  preceding  discussion  assumed  that  the  theoret¬ 
ical  curves  are  accurate.  Is  this  assumption  valid?  The 
droplet  trajectory  analysis  upon  which  the  curves  are 
based  was  done  originally  for  an  incompressible  flow 
field  with  the  differential  analyzer  at  the  General  Elec¬ 
tric  Research  Laboratory  (Langmuirand  Blodgett  1 946). 
This  analysis  was  redone  for  a  compressible  flow  field 
by  NACA,  using  a  mechanical  analog  based  on  the  di¬ 
fferential  analyzer  (Brun  et  al.  1955).  The  NACA  re¬ 
sults  were  slightly  different,  and  NACA  also  improved 
the  method  for  computing  the  RMC  curves.  We  use  the 
NACA  data  and  computation  method.  The  trajectory 
analysis  was  repeated  on  a  computer  by  McComber  and 
Touzot  (1981)  using  the  method  of  finite  elements. 
Their  work  confirmed  the  accuracy  of  the  NACA  re¬ 
sults;  there  were  some  differences,  but  these  would  not 
alter  the  RMC  curves  significantly.  Although  some  un¬ 


certainty  may  remain  regarding  the  effects  of  turbu¬ 
lence  and  of  cylinder  surface  roughness,  there  is  a  high 
degree  of  confidence  in  the  accuracy  of  the  RMC  curves 
for  monodisperse  droplet  populations. 

A  fundamental  and  potentially  serious  problem  re¬ 
mains,  however,  concerning  the  development  of  the 
curves  fbr  the  B  through  J  droplet  size  distributions. 
Langmuir’s  original  method  has,  to  my  knowledge, 
been  used  in  all  subsequent  work.  (My  own  attempts  at 
revision  will  be  described  below.)  Langmuir  adopted 
one  of  Houghton  and  Radford’ s  (1938)  measured  distri¬ 
butions  of  fog  droplet  size  as  a  reasonable  basis  for  what 
he  defined  as  the  B  distribution;  his  distributions  C,  D, 
and  E  were  derived  mathematically  therefrom,  and  the 
broader  distributions  were  added  by  the  MWO  using  the 
same  procedure.  These  distributions  are  approximately 
normal  (Gaussian),  and  can  be  represented  as  histo¬ 
grams  of  total  droplet  volume  vs  droplet  diameter  for 
equal  increments  of  droplet  size. 

It  was  recognized  during  the  first  seasons  of  routine 
RMC  exposures  at  the  MWO  that  bimodal  distributions 
must  occasionally  occur,  particularly  during  episodes 
of  fog  with  freezing  drizzle  or  rain.  Theoretical  curves 
representing  some  hypothetical  bimodal  distributions 
were  produced  (these  generally  recurve,  becoming  con¬ 
cave  upwards  for  the  higher  values  of  1/AT),  but  these 
curves  were  not  found  to  be  very  useful.  No  doubt  this 
attempt  failed  simply  because  the  possible  variations  in 
nature  of  such  bimodal  distributions  are  practically 
infinite. 

I  approached  this  problem  in  the  early  1960s  by 
synthesizing  sets  of  data  points  for  RMC  exposures  in 
hypothetical  bimodal  and  trimodal  clouds.  In  all  but  one 
of  the  seven  cases  studied,  the  synthesized  data  points 
could  be  fitted  to  one  of  the  standard  theoretical  curves. 
Table  1  summarizes  these  results.  Note  that  in  all  but  the 
third  case  the  resulting  value  of  LWC  is  within  10%  of 
the  sum  of  the  LWCs  in  the  synthesized  populations, 
and  the  resulting  volume  median  droplet  diameter  and 
distribution  are  at  least  reasonable. 

During  the  1970s  the  MWO  made  a  series  of  simul¬ 
taneous  exposures  with  the  RMC,  the  Blau  Optical 
Cloud  Particle  Spectrometer(OCPS)  (Ryan  et  al.  1 972), 
and  an  oiled-slide  droplet  collector.  (Results  of  these 
tests  are  summarized  in  a  later  section.)  The  raw  data 
from  the  latter  two  instruments  were  in  the  form  of 
droplet  counts,  and  histograms  in  most  ases  showed  an 
approximately  normal  or  Gaussian  distr:  bution  c  f  drop¬ 
let  count  over  the  range  of  droplet  size.  This  result  was 
unexpected;  whether  or  not  it  was  entirely  correct,  it  at 
least  strongly  suggested  that  the  Houghton/Langmuir 
distributions,  approximately  normal  with  respact  to 
droplet  volume,  might  not  be  representative  of  condi¬ 
tions  in  the  real  world. 
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Table  1.  Synthesized  RMC  data. 


LIVC 

Droplet  populations. 

synthesized  data 

Results  using 
standard  procedure 
LWC  Dd  Dist. 

/ 

D/ 

Disl. 

LWC 

II 

D, 

Disl. 

LWC 

Ill 

Dist. 

0.1 

8 

A 

0.2 

40 

A 

_ 

_ 

_ 

0.27 

32 

F 

0.1 

8 

D 

0.1 

40 

D 

— 

— 

— 

0.19 

19 

J 

0.1 

8 

D 

0.3 

100 

D 

— 

— 

— 

0.31 

90 

E 

0.1 

8 

A 

0.1 

20 

A 

0.2 

40 

A 

0.38 

25 

E 

0.1 

8 

D 

0.03 

40 

D 

— 

— 

— 

0.15 

9 

J 

0.1 

8 

D 

0.03 

100 

D 

— 

— 

— 

(0.12) 

Poor  Fit 

0.2 

20 

B 

0.03 

100 

B 

— 

— 

— 

0.23 

23 

D 

I  therefore  constructed  new  sets  of  theoretical  curves 
with  a  distribution  scheme  based  on  a  typical  OCPS 
exposure.  For  the  past  1 4  years  most  RMC  runs  made  at 
the  MWO  (several  hundred  runs)  have  been  analyzed  on 
both  the  new  and  the  old  (NACA/Langmuir)  curves.  In 
every  case,  there  was  very  little  difference  in  goodness 
of  fit.  Results  generally  differed  by  less  than  5%  in 
LWC  and  Dd  for  narrow  and  medium  distributions:  with 
broad  distributions  differences  were  usually  less  than 
10%,  but  occasionally  as  much  as  20%.  The  new  curves 
usually,  but  not  always,  gave  higher  values.  Indicated 
distributions  were  always  broader  with  the  new  curves, 
as  would  be  expected. 

The  two  investigations  described  above  show  that 
the  RMC  method  is  insensitive  to  the  details  of  droplet 
size  distribution.  The  fact  that  a  data  plot  fits  a  particular 
theoretical  curve  does  not  imply  that  the  cloud  had  the 
size  distribution  assumed  for  the  curve,  and  it  may  in 
fact  have  been  quite  different.  The  method  gives  only  a 
general  indication  of  the  breadth  of  droplet  size  distri¬ 
bution.  Fortunately,  however,  this  shortcoming  has  a 
rather  small  effect  on  the  accuracy  of  LWC  and  Dd 
values  except  when  the  distribution  is  very'  broad. 

Comparison  tests 

Howell  (1952)  investigated  the  performance  of  three 
RMC  instruments  that  differed  in  several  design  details, 
and  found  no  significant  differences  in  results.  I  did  a 
few  such  tests  with  two  instruments  in  1985.  with  the 
same  conclusion. 

The  simultaneous  exposures  of  RMC.  OCPS.  and  an 
oiled-slide  collector,  mentioned  above,  show  ed  reason¬ 
ably  good  correlation  in  £>d,  with  OCPS  values  ap¬ 
proximately  20%  higher  than  RMC  values.  This  may  be 
partly  explained  as  being  due  to  the  presence  during 
most  runs  of  snow  or  blowing  snow .  which  affected  the 
OCPS  but  not  the  RMC.  The  OCPS  gave  LWCs  that 
were  usually  half  to  one-tenth  of  those  given  by  the 
RMC;  no  explanation  for  this  was  found  despite  much 
effort.  No  attempt  was  made  to  compute  LWC  from  the 
oiled-slide  data. 


Simultaneous  exposures  of  the  RMC  and  the  Parti¬ 
cle  Measuring  Systems  Forward  Scattering  Spectrom¬ 
eter  Probe  (FSSP)  at  the  MWO  during  1984-85  gave 
similar  results,  i.e.,  fairly  good  agreement  inDd,  but  FSSP 
values  of  LWC  that  were  much  too  small. 

The  large  discrepancies  in  LWC  values  are  certainly 
due  to  poorly  understood  errors  in  the  optical  probes' 
sensing  and/or  data  processing.  This  statement  sounds 
arbitrary,  but  consider  the  following: 

1.  RMC  physical  measurements  are  straightforward 
and  unlikely  to  be  in  error  by  more  than  10%. 

2.  Even  if  all  the  arguments  given  in  preceding 
sections  are  disregarded  for  accuracy  of  RMC  data 
analysis  methods,  the  EW  value  measured  for  the 
smallest  cylinder  represents  a  lower  limit  for  the  value 
of  LWC.  The  smallest-cylinderEW  values  forthese  runs 
were  generally  60%  to  80%  of  the  LWC  values. 

It  should  be  mentioned  that  both  optical  probes  were 
calibrated  immediately  prior  to  the  MWO  tests.  The 
calibration  of  the  OCPS  was  checked  immediately  after 
the  tests  and  found  to  be  unchanged. 


CONCLUSIONS 

The  Rotating  Multicylinder  (RMC)  is  a  low-cost, 
simple,  reliable  instrument.  Exposure  and  data-reduction 
procedures  are  straightforward;  they  are  somewhat  time- 
consuming,  but  the  almost  total  lack  of  maintenance  and 
repair  requirements  largely  offsets  this  disadvantage. 
Although  intended  for  use  in  supercooled  clouds,  the 
method  can,  with  some  modifications,  be  used  in  above¬ 
freezing  temperatures.  With  reasonable  care  in  making 
the  necessary  physical  measurements,  accuracy  in  de¬ 
termination  of  liquid-water  content  and  droplet  size  is 
better  than  ±10%  when  cloud  droplet-size  distribution 
is  narrow  or  moderately  broad.  For  extremely  broad 
distributions,  accuracy  may  be  reduced  to  about  ±20%. 
The  method  is  insensitive  to  details  of  size  distribution, 
but  it  reliably  indicates  the  approximate  distribution 
breadth. 
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This  accuracy  is  achieved  partly  because  the  method 
samples  an  astronomically  large  number  of  droplets 
over  a  large  volume  of  cloud.  This  may  be  a  disadvan¬ 
tage  in  some  applications. 

Errors  in  measurement  will  generally  be  revealed 
during  data  reduction  and  can  easily  be  allowed  for;  in 
most  such  cases  accuracy  will  not  be  reduced.  Accuracy 
is  worst  when  the  temperature  is  a  few  degrees  below 
freezing,  due  to  difficulties  in  handling  and  the  likeli¬ 
hood  of  runoff. 
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APPENDIX  A:  MULTICYLINDER  INSTRUMENT  DESIGN 


The  shaft  of  the  RMC  should  be  made  of  stainless 
steel  or  brass  and  the  cylinders  and  intercylinder  adap¬ 
tors  of  aluminum,  except  for  the  two  smallest  cylinders, 
which  may  more  easily  be  made  of  brass.  The  diameter 
of  the  cylinders  is  not  critical:  standard  tubing  sizes  can 
be  used  such  that  each  cylinder  is  approximately  twice 
the  diameter  of  its  smaller  neighbor  and  the  diameter  of 
the  largest  is  about  3  in.  (7.62  cm). 

A  slightly  loose  fit  between  all  the  components 
makes  for  ease  of  disassembly.  The  bottom  adaptor  is 
secured  to  the  shaft  with  epoxy.  The  no.  2  cylinder  has 
a  flange  and  threaded  lower  end,  and  screws  into  the  top 
end  of  the  shaft  to  hold  every  thing  i  n  place.  The  smallest 
cylinder  fits  into  a  drilled  hole  in  the  top  of  the  no.  2 
cylinder. 

Figure  Al  shows  the  design  for  the  shaft  and  adap¬ 
tors.  Dimensions  are  given  in  English  units  for  simplic¬ 
ity:  in  the  U.S.  it  will  be  easiest  to  obtain  the  tubing  for 
the  cylinders  in  standard  inch  dimensions.  Note  that 
dimensions  A  and  F  for  the  adaptors  are  equivalent  to 


the  associated  cylinder  diameters,  and  that  B  and  E  will 
depend  on  the  wall  thickness  of  the  tubing  used  for  the 
cylinders.  The  flanges  shown  on  the  adaptors  are  op¬ 
tional;  our  experience  is  that  they  do  no  harm  and  are 
helpful  at  times  during  disassembly.  They  influence  the 
ice  collection  only  if  the  RMC  axis  is  not  perpendicular 
to  the  airstream,  and  in  fact  provide  a  very  sensitive 
indication  of  that  undesirable  condition. 

At  the  MWO  our  no.  2  cylinder  is  machined  from 
solid  brass  rod.  Finished  length  is  1 1  cm,  diameter  -0.5 
cm.  The  bottom  1  cm  is  threaded  to  fit  the  tapped  hole 
in  the  top  of  the  shaft  with  a  small  flange  -0.6  cm  in 
diameter  to  bear  against  the  top  of  the  1I/II1  adaptor.  Our 
no.  1  cylinder  is  just  a  piece  of  0.158-cm  brass  wire 
about  10  cm  long;  the  bottom  end  can  be  crimped 
slightly  so  that  it  fits  easily  but  securely  in  the  drilled 
hole  at  the  top  end  of  the  no.  2  cylinder.  The  four  largest 
cylinders  are  made  from  aluminum  tubing  with  ends 
machined  to  give  a  length  of  7.4  cm. 


II 


APPENDIX  B:  THEORETICAL  CURVES 


Droplet  size  distributions 

The  droplet  size  distributions  assumed  for  the  theo¬ 
retical  curves  are  listed  in  Table  B 1 .  These  are  the  same 
distributions  assumed  in  Langmuir  and  Blodgett  (1946) 
and  in  Brun  et  al.  (1955),  except  for  the  addition  of  the 
broader  types,  F,  G,  H,  and  J. 

In  Table  B 1  the  droplets  in  each  distribution  type  are 
divided  into  seven  groups  representing  symmetrical 
subdivisions  of  the  area  under  the  distribution  curve. 
The  values  ala d  are  the  ratios  of  the  droplet  radius  for 
each  size  group  to  the  radius  of  the  volume-median  size. 
The  volume  fraction  of  liquid  water  in  each  size  group 
is  listed  in  the  left  column.  For  the  A  distribution  all  the 
droplets  are  the  same  size.  The  values  of  a!a&  for  the  B 
distribution  are  raised  to  the  1 .5,  2.0,  2.5,  3.0,  3.5, 4.0, 
and  5.0  power  to  obtain  the  ala ^  values  for  the  C,  D,  E, 
F,  G,  H,  and  J  distributions,  respectively. 

Calculation  of  curves  of  UK  vs  Em 

The  theoretical  curves  of  1  IK  vs  Em  shown  in  Figure 
B1  are  constructed  from  the  values  given  in  Table  B3. 
The  method  of  Brun  et  al.  (1955)  was  used  to  calculate 


these  values;  they  gave  tables  of  values  for  K§  =  0, 200, 
1000, 3000,  and  10,000  and  distributions  A  through  E, 
which  are  reproduced  here  as  part  of  Table  B3.  The  re¬ 
maining  values  were  calculated  by  me  at  the  Aeronaut¬ 
ical  Icing  Research  Laboratories  in  1 956.  At  that  time  I 
also  calculated  values  for  1  IK  =  10.0  for  all  Kt\>  values 
and  distributions  where  necessary  to  carry  the  curves  to 
the  top  edges  of  the  graphs.  (This  was  not  done  by  Brun 
et  al.). 

Values  for  each  of  the  A  curves  are  calculated  by 
assigning  a  series  of  values  to  \/K  and  finding  the  cor¬ 
responding  values  of  £  from  the  curves  of  E  vs  K,  with 
(j)  as  a  parameter  given  by  Brun  et  al.  (The  value  of  <J>  is 
determined  for  each  value  of  l/K,  since  £0  is  held 
constant  for  each  curve  of  UK  vs  Em-). 

For  the  B  through  J  distributions  of  Table  B3.  overall 
weighted  collection  efficiencies  (Em)  aie  used.  Table 
B2,  which  shows  the  procedure  for  finding  a  single 
point  of  the  K§  =  100-  J  distribution  curve,  will  serve  as 
an  example .  Here  the  assigned  value  of  1  /K  is  0.0 1  (K  for 
volume-median  size  equals  1 00. 0  =  1 0).  For  each  size 
group,  the  appropriate  value  of  K  is  found  by  multiply- 
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Table  BI.  Droplet  size  distributions. 


Fraction  of 

LWC  in  each  Distribution  (a/ad) 


size  group 

A 

B 

C 

D 

£ 

F 

G 

H 

1 

0.05 

1.00 

0.56 

0.42 

0.31 

0.23 

0.18 

0.13 

0.10 

0.56 

0.10 

1.00 

0.72 

0.61 

0.52 

0.44 

0.37 

0.32 

0.27 

0.194 

0.20 

1.00 

0.84 

0.77 

0.71 

0.65 

0.59 

0.54 

0.50 

0.42 

0.20 

1.00 

1.17 

1.26 

1.37 

1.48 

1.60 

1.73 

1.88 

2.20 

0.30 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

o 

p 

0.10 

1.00 

1.32 

1.51 

1.74 

2.00 

2.30 

2.64 

3.03 

4.00 

0.05 

1.00 

1.49 

1.81 

2.22 

2.71 

3.31 

4.04 

4.93 

7.34 

ing  the  vol  jme-median  size  lvalue  by  the  value  of  (a/ 
ad)2  for  that  size  group,  since  the  droplet  size  is  raised 
to  the  second  power  in  the  expression  for#.  The  value 
of  <j>  remains  the  same  as  for  the  volume-median  size 
since  droplet  size  does  not  appear  in  the  expression  for 
ij).  The  collection  efficiency  for  each  size  group  is  then 
obtained  from  the  curves  of  Brun  et  al.  (1955)  men¬ 
tioned  above.  The  weighted  collection  efficiency  for 
each  size  group  is  obtained  by  multiplying  the  value  in 
column  three  by  the  value  in  column  five.  The  sum  of 
column  six  is  the  overall  weighted  collection  efficiency, 
givingthe  point  1/A" = 0.0 1  ,£TM =0.899  for  the  /f<j)  =  1000 
J  curve. 

This  procedure  differs  from  that  of  Langmuir  and 
Blodgett  (1946)  in  that  the  value  of  <|>  is  held  constant 
foreach  point,  whereas  in  Langmuir  and  Blodgett  <{>  was 
varied  to  keep  K<f>  constant.  The  method  of  Brun  et  al.  is 
correct,  since  <i>  is  not  a  function  of  D. 

Figure  B 1  shows  an  example  of  the  curves  used  at 
MWO.  All  but  a  tiny  fraction  of  exposures  in  natural 


wind  will  require  only  the  curves  for  K<f>  =  20, 200,  and 
1000. 

Curve  matching 

Brun  et  al.  (1955)  give  a  rigorous  exposition  of  the 
procedure  for  curve-matching  and  the  determination  of 
W  and  CK_  j .  Briefly,  the  procedure  can  be  explained  by 
stating  that  the  multiplication  and  division  of  numbers 
required  for  these  determinations  is  equivalent  to  the 
addition  and  subtraction  of  the  logarithms  of  those 
numbers,  and  this  is  done  graphically  by  sliding  the  log- 
log  graph  of  plotted  points  parallel  to  the  theoretical 
curves,  which  are  on  the  same  type  of  log-log  paper. 

For  a  more  intuitive  understanding,  note  that  the 
theoretical  curves  are  asymptotic  to  the  £M  =  1  axis  for 
small  valuesof  l/K;  the  inverted  T  that  you  traceon  your 
data  plot  is  the  location  of  that  axis.  Stating  it  another 
way,  if  your  smallest  cylinder  had  a  collection  efficien¬ 
cy  of  100%  (£=  1),  its  data  point  would  have  fallen  on 
that  axis  and  its  EW  value  would  have  been  equal  to  W. 


Table  B2.  Sample  calculation  of  EM:  K§  =  1000,  K  =  100, 
distribution  J. 


/ 

2 

3 

4 

5 

6 

Fraction  of 

Weighted 

Sire 

water  in 

collection 

group 

(alad)2 

size  group 

K 

E 

efficiency 

i 

0.003 1 

0.05 

0.31 

0.087 

0.0044 

2 

0.0377 

0.10 

3.77 

0.715 

0.0715 

3 

0.177 

0.20 

17.7 

0.9)4 

0.1828 

4 

1.00 

0.30 

100.0 

0.972 

0.2916 

5 

4.84 

0.20 

484.0 

0.995 

0.1990 

6 

16.00 

0.10 

1600.0 

0.999 

0.0999 

7 

54.0 

0.05 

5400.0 

0.999 

0.0500 

Overall  weighted  collection  efficiency.  £M  =  0.899 
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Table  B3.  Data  for  construction  of  theoretical  curves. 


K6  IK 

A 

B 

C 

D 

£ 

F 

G 

H 

/ 

0  10 

_ 

_ 

0.010 

0.022 

0.036 

0.053 

0.072 

0.091 

0.130 

4 

0.051 

0.069 

0.085 

0.107 

0.126 

0.146 

0.167 

0.186 

0.221 

2 

0.205 

0.204 

0.211 

0.225 

0.241 

0.256 

0.273 

0.290 

0.318 

i 

0.380 

0.374 

0.373 

0.379 

0.384 

0.389 

0.395 

0.402 

0.413 

0.5 

0.566 

0.555 

0.549 

0.542 

0.536 

0.529 

0.523 

0.520 

0.514 

0.2 

0.789 

0.768 

0.750 

0.732 

0.713 

0.695 

0.679 

0.668 

0.643 

0.1 

0.885 

0.870 

0.854 

0.836 

0.815 

0.794 

0.773 

0.758 

0.727 

0.05 

0.932 

0.925 

0.918 

0.898 

0.885 

0.867 

0.846 

0.828 

0.798 

0.02 

0.963 

0.961 

0.957 

0.951 

0.940 

0.929 

0.913 

0.899 

0.866 

0.01 

0.978 

0.976 

0.977 

0.972 

0.965 

0.958 

0  946 

0.934 

0.905 

20  10 

_ 

0.0056 

0.0137 

0.0247 

0.038 

0.056 

0.069 

0.104 

4 

0.038 

0.052 

0.067 

0.085 

0.103 

0.122 

0.142 

0.161 

0.197 

2 

0.170 

0.173 

0.181 

0.195 

0.212 

0.228 

0.246 

0.262 

0.293 

i 

0.342 

0.340 

0.341 

0.346 

0.354 

0.362 

0.370 

0.376 

0.389 

0.5 

0.540 

0.537 

0.524 

0.520 

0.514 

0.510 

0.505 

0.503 

0.499 

0.2 

0.780 

0.759 

0.739 

0.721 

0.698 

0.682 

0.667 

0.655 

0.632 

0.1 

0.880 

0.864 

0.849 

0.830 

0.806 

0.787 

0.767 

0.752 

0.720 

0.05 

0.930 

0.921 

0.913 

0.898 

0.882 

0.865 

0.842 

0.827 

0.796 

0.02 

0.962 

0.959 

0.956 

0.949 

0.940 

0.928 

0.913 

0.898 

0.865 

0.01 

0.978 

0.976 

0.974 

0.971 

0.966 

0.957 

0.946 

0.934 

0.905 

200  10 

— 

0.0037 

0.0092 

0.0168 

0.0269 

0.040 

0.054 

0.086 

4 

0.027 

0.039 

0.050 

0.066 

0.083 

0.100 

0.118 

0.137 

0.173 

2 

0.135 

0.138 

0.146 

0.165 

0.182 

0.196 

0.216 

0.229 

0.261 

1 

0.298 

0.302 

0.306 

0.315 

0.319 

0.323 

0.334 

0.345 

0.361 

0.5 

0.493 

0.486 

0.482 

0.480 

0.477 

0.474 

0.471 

0.472 

0.471 

0.2 

0.761 

0.740 

0.721 

0.703 

0.686 

0.665 

0.646 

0.635 

0.615 

0.1 

0.874 

0.859 

0.846 

0.826 

0.805 

0.777 

0.754 

0.739 

0.705 

0.05 

0.925 

0.919 

0.910 

0.901 

0.878 

0.862 

0.841 

0.819 

0.783 

0.02 

0.960 

0.959 

0.955 

0.948 

0.938 

0.927 

0.911 

0.894 

0.866 

0.976 

lilCMV 

■ifcjcn 

■ivjixa 

10-’  10 

— 

— 

0.0066 

0.0121 

0.0193 

0.028 

0.039 

0.066 

4 

0.019 

0.029 

0.038 

0.050 

0.065 

0.076 

0.092 

0.110 

0.142 

2 

0.109 

0.109 

0.122 

0.138 

0.154 

0.164 

0.181 

0.199 

0.231 

i 

0.251 

0.252 

0.259 

0.271 

0.283 

0.293 

0.304 

0.315 

0.334 

0.5 

0.460 

0.452 

0.423 

0.447 

0.447 

0.438 

0.437 

0.438 

0.440 

0.2 

0.714 

0.697 

0.677 

0.661 

0.643 

0.623 

0.609 

0.601 

0.585 

0.1 

0.830 

0.816 

0.800 

0.783 

0.763 

0.743 

0.724 

0.711 

0.683 

0.05 

0.908 

0.899 

0.892 

0.876 

0.862 

0.835 

0.815 

0.800 

0.768 

0.02 

0.953 

0.953 

0.949 

0.943 

0.933 

0.912 

0.896 

0.882 

0.851 

0  01 

0  971 

n  976 

0  973 

0  967 

0  937 

0  977 

5x1  O'*  10 

0.0041 

0.0088 

0.014 

0.021 

0.031 

0.054 

4 

0.013 

0.020 

0.027 

0.039 

0.048 

0.060 

0.075 

0.088 

0.120 

2 

0.085 

0.090 

0.100 

0.111 

0.130 

0.140 

0.159 

0.177 

0.211 

1 

0.218 

0.225 

0.235 

0.244 

0.251 

0.256 

0.268 

0.282 

0.303 

0.5 

0.409 

0.416 

0.410 

0.415 

0.413 

0.408 

0.399 

0.405 

0.411 

0.2 

0.687 

0.668 

0.652 

0.641 

0.623 

0.601 

0.591 

0.583 

0.567 

0.1 

0.815 

0.797 

0.785 

0.766 

0.746 

0.719 

0.702 

0.689 

0.663 

0.05 

0.884 

0.878 

0.867 

0.855 

0.839 

0.813 

0.795 

0.779 

0.750 

0.02 

0.954 

0.940 

0.938 

0.921 

0.918 

0.901 

0.885 

0.870 

0.840 

nm 

0  968 

0  966 

0  970 

0964  .. 

0.954 

0  940 

0.977 

0917 

0.8S8 

104  10 

0.0082 

0.0136 

0.0207 

0.040 

4 

0.008 

0.013 

0.017 

0.023 

0.034 

0.046 

0.059 

0.073 

0.103 

2 

0.057 

0.060 

0.072 

0.083 

0.092 

0.110 

0.126 

0.144 

0.176 

i 

0.157 

0.163 

0.172 

0.188 

0.202 

0.216 

0.231 

0.246 

0.272 

0.5 

0.350 

0.356 

0.357 

0.362 

0.368 

0.356 

0.366 

0.374 

0.385 

0.2 

0.645 

0.630 

0.615 

0.599 

0.591 

0.559 

0.550 

0.542 

0.533 

0.1 

0.778 

0.764 

0.748 

0.731 

0.713 

0.695 

0.677 

0.666 

0.643 

0.05 

0.865 

0.857 

0.849 

0.830 

0.816 

0.793 

0.777 

0.762 

0.734 

0.02 

0.920 

0.920 

0.918 

0.909 

0.899 

0.884 

0.869 

0.855 

0.825 

..  0.0 L 

_ 0.952 

_ 0.950 _ 

0.955 

0.946 

0.939 

0.930 

0.918 

0.907 

0.879 

UP  10 

4 

— 

0.0045 

0.0073 

0.0110 

0.00156 

0.0202 

0.0307 

0.0391 

0.061 

2 

0.022 

0.026 

0.031 

0.038 

0.048 

0.060 

0.073 

0.089 

0.122 

1 

0.080 

0.086 

0.094 

0.109 

0.125 

0.140 

0.156 

0.171 

0.200 

0.5 

0.225 

0.226 

0.233 

0.243 

0.254 

0.266 

0.277 

0.290 

0.310 

0.2 

0.485 

0.479 

0.474 

0.470 

0.465 

0.463 

0.462 

0.463 

0.460 

0.1 

0.687 

0.669 

0.654 

0.640 

0.625 

0.612 

0.599 

0.593 

0.578 

0.5 

0.805 

0.794 

0.782 

0.768 

0.750 

0.734 

0.719 

0.706 

0.681 

0.02 

0.887 

0.883 

0.875 

0.870 

0.858 

0.846 

0.830 

0.818 

0.789 

0.01 

0.930 

0.926 

_ 0524 

0.920 

0.913 _ 

0.904 

Q.891 _ 

0.880 

0.855 
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APPENDIX  C:  STEP-BY-STEP  PROCEDURE 


The  numbers  in  Table  C 1  are  taken  from  an  exposure 
made  at  MWO  at  0332  on  7  Dec  1987.  The  value  of  Wj 
for  a  no.  3  cylinder  has  been  changed  to  simulate  an 
error  in  weighing.  Some  details  of  the  procedure  have 
been  omitted,  since  they  will  depend  on  individual 
circumstances  (exposure  site,  computer  usage,  etc.). 
One  important  service  of  this  Appendix  is  to  demon¬ 
strate  the  magnitude  of  the  numbers  you  can  expect  to 
measure  and  compute;  without  such  help  it  may  be 
difficult  to  keep  track  of  the  decimal  point  since  the  unit 
of  length  is  variously  meters,  centimeters,  or  microme¬ 
ters. 

Exposure 

1 .  Make  sure  the  instrument  is  dry  and  cold.  If  it  has 
been  cooled  while  still  wet  it  may  be  impossible  to 
disassemble  it  after  the  exposure,  so  check  to  see  if  all 
the  parts  move  freely.  Also  check  the  exposure  site  and 
prepare  to  deploy  your  anemometer  or  deice  a  perma¬ 
nently  mounted  one. 

2.  Position  the  instrument  for  exposure  and  punch 
your  stopwatch.  Then,  for  a  few  moments,  observe 
what's  going  on.  Is  the  RMC  rotating?  If  conditions  are 
marginal,  is  water  running  down  the  smaller  cylinders 

Table  Cl.  Sample  exposure  data  and  analysis. 

Pressure  (mb):  809  Temp  (°C):  -14  Wind  speed  (mph):  66 
Exposure  duration(s):  900  pa  (g/cm-'):  1.075  x  10~3 

V/r\-  18.02  x  106  /,:  _187 


Cylinder 


no. 

1 

2 

3 

4 

5 

6 

C, 

0.462 

0.758 

1.33 

_ 

m 

— 

— 

— 

0.586 

0.560 

0.568 

Co 

0.158 

0.502 

1.11 

— 

— 

— 

n 

o 

1 

— 

— 

— 

1.738 

3.955 

7.115 

0.310 

0.630 

1.22 

2.32 

4.52 

7.68 

L 

9.20 

10.25 

7.35 

7.40 

7.35 

7.40 

wi 

2.84 

19.46 

8.08 

21.31 

42.13 

78.88 

“V 

1.84 

17.64 

5.79 

18.90 

39.72 

77.03 

EW 

0.132 

0.106 

0.0962 

0.0529 

0.0273 

0.0123 

CK=,  (cm):  Wlg/m-1):  017 

K< t>:  542  0d(gm>:  11 

Curve  used:  tOOO  Distribution:  (7 

Points  fitted:  5 


without  freezing?  (If  so,  give  up  and  try  again  another 
day.)  Check  the  progress  of  the  run  often  enough  to 
avoid  excessive  accumulation. 

3.  End  the  exposure  when  the  deposit  on  the  smallest 
cylinder  has  grown  to  approximately  the  original  diam¬ 
eter  of  the  next-larger  one.  Timing  to  within  about  five 
seconds  is  adequate.  Bring  the  instrument  back  into  the 
coldroom.  Be  careful;  sooner  or  later  you  will  probably 
be  careless  and  break  some  ice  or  even  one  of  the  two 
smallest  cylinders.  (You  should  have  spares  for  those 
small  cylinders.)  Place  the  instrument  in  a  prepared 
position  so  that  it  stands  securely  upright.  With  your 
micrometer  caliper,  measure  the  diameters  of  the  three 
smallest  cylinders.  If  the  collection  is  fragile  rime,  be 
careful  not  to  crush  it;  adjust  the  micrometer  by  eye,  not 
by  feel.  Remember  to  enter  the  duration  lime  and  all 
other  measurements  on  your  work  sheet,  and  also  indi¬ 
cate  the  characterof  the  accumulation  on  each  cylinder. 

4.  Put  on  apairof  cotton  gloves,  remove  the  smallest 
cylinder,  measure  the  length  of  the  collection  to  the 
nearest  0.05  cm,  and  weigh  to  the  nearest  0.0 1  g.  Repeat 
for  the  other  cylinders,  one  at  a  time.  If  the  ice  is  be¬ 
ginning  to  melt,  weigh  each  cylinder  before  you  meas¬ 
ure  its  diameter.  Diameters  of  the  three  largest  ones  are 
measured  indirectly  over  a  wire  anvil  (see  last  para¬ 
graph  under  Instrumentation).  Separating  the  cylinders 
from  the  intercylinder  adapters  will  be  difficult  if  the 
collection  is  clear  ice  and  you  will  have  to  develop  your 
own  technique,  but  see  under  Exposure  Procedure  for 
a  suggestion. 

Data  reduction 

5.  Compute  the  average  diameter  Cav  for  each  cyl¬ 
inder.  For  the  three  smallest  ones  this  is  simply  equal  to 
(Cj  +  C0)/ 2.  For  the  three  largest  cylinders,  which  were 
measured  over  an  anvil,  use  the  following  expression: 

Cav  =  m  +  ( C-b ) 

where  m  -  (anvil)+(cylinder  wall)+(ice)  and  b  = 
(anvil)+(cylinder  wall). 

6.  Compute  EW  for  each  cylinder  by  eq  1  in  body  of 
this  report  (Y  in  cm/s,  L  and  Cav  in  cm). 

7.  Compute  Wr|  and  _/j  by  eq  3  and  4.  Note  that 
temperature  must  be  in  kelvins  (=  °C+273,  e.g..  -14°C 
=  295  K). 

8.  On  three-cycle  log  paper,  plot  the  values  of  Cav 
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and  EW,  with  Cav  on  the  vertical  axis  and  EW  on  the 
horizontal  axis  (Fig.  2a).  (You  will  have  to  choose  the 
values  for  the  axes  such  that  the  plot  will  be  near  the 
center  of  the  paper.) 

9.  Now  for  the  curve-matching.  You  have  to  guess 
which  family  of  theoretical  curves  to  try  first  (usually 
AT4>  =  20,  200.  or  1 000).  Say,  for  this  example,  that  you 
pick  the  200  family.  Keeping  the  axes  parallel,  slide 
your  plotted  graph  around  on  the  theoretical  curves  until 
you  find  the  best  fit  (in  this  case  the  C  curve).  You  will 
see  an  inverted  T  on  the  theoretical  curve,  labeled  C, 
with  its  horizontal  arm  on  the  \/K=  1  axis.  Note  that  this 
is  under  the  value  Cav  =  2.5  on  your  plot,  meaning  that 
CK=1  =2.5.  Multiply  this  value  by  the  value  of /j  which 
you  found  in  Step  7  (2.5x187  =  468).  Determine 
whether  468  is  closer  to  200  than  it  is  to  1000,  on  a 
logarithmic  scale  (you  can  use  the  scale  on  the  graph 
paper  for  this).  The  answer  is  “no,”  it’s  closer  to  1000. 
This  means  that  you  must  repeat  the  matching  on  the  K(j> 
=  1000  curves. 

1 0.  Do  your  final  curve-matching  on  the  1 000  curves 


and  you  will  find  that  the  best  fit  is  again  on  the  C  curve. 
Lightly  trace  the  curve  through  your  plotted  points  and 
also  trace  the  inverted  T  labeled  C.  (Its  location  will  be 
slightly  different  from  what  it  was  when  you  tried  the 
200  curve.)  Read  the  value  of  the  horizontal  arm  of  the 
T  on  your  graph  (CK_,  =2.9).  Then  read  the  value  of  the 
vertical  arm  of  the  Ton  your  graph  (£WE=1  =0.1 7).  This 
means  that  the  liquid-water  content  (LWC)  for  your 
exposure  was  0.17  g/m. 

1 1 .  Multiply  your  value  of  CK= ,  by ,/j  to  find  the  value 
of  AT0.  Calculate  the  value  of  Dd  by  eq  5.  In  the  space  for 
“Distribution”  on  your  analysis  sheet,  enter  C,  and 
indicate  the  number  of  data  points  which  fit  the  theoret¬ 
ical  curve  within  5%. 

Important  note:  It  is  essential  that  the  theoretical 
curves  be  on  the  same  type  of  graph  paper  which  you  use 
for  your  plots  of  Cav  vs  EW.  You  must  construct  your 
own  curves  from  the  values  given  in  Table  B3.  It  is 
important  to  use  three-cycle  by  three-cycle  paper  with 
as  large  a  format  as  practicable.  Construction  requires 
skillful  use  of  french  curves. 
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